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ABSTRACT 



A five degree-of-freedom linear model for the XR-3 
surface effect ship is developed for constant speed operation. 
The weight removal transient response and the roll behavior 
of the craft in turn maneuvers are investigated by both 
linear and six degree-of-freedom nonlinear models, and the 
linear model simulation results are compared with the non- 
linear model simulations. Some of the parameters of the 
craft are investigated for roll motion. 
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I . INTRODUCTION 



A. BACKGROUND 

The Surface Effect hybrid type high performance ship has 
been of great interest recently, and in 1970 the XR-3 test 
craft was delivered to the Naval Postgraduate School. As 
seen in [7] Leo and Boncal converted the digital simulation 
program for the 100-B test craft to represent the XR-3, thus 
providing a nonlinear model. In [1] , Gerba and Thaler have 
developed a heave-only linear model and further studies were 
made by Barnes in [2] by developing a two degree-of- freedom 
linear pitch and heave model representing the XR-3 test craft. 

B. OBJECTIVE 

In order to investigate the roll behavior of the XR-3 
test craft during turn maneuvers, a five degree-of-freedom 
linear model is developed in this thesis. Some modifications 
are made for pitch and heave modes, therefore test of the 
vertical plane by means of weight removal is reinvestigated 
and compared with the nonlinear simulation results. 

In [3] , the roll behavior of the XR-3 and the sidewall 
sensitivity are studied using a six degree-of-freedom non- 
linear model. Similar studies are made in this thesis after 
developing a five degree-of-freedom linear model and comparing 
with the nonlinear model results. Finally, some modifications 
are suggested for nonlinear modeling. 
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II. NONLINEAR XR-3 MODELING FOR FIVE DEGREE-OF-FREEDOM 



A. ASSUMPTIONS 

The following assumptions and simplifications are used 
during development of a constant speed XR-3 test craft model; 

1. Aerodynamic drag on the forward seal, and superstruc- 
ture skin friction forces are neglected. 

2. Momentxim (ram) drag is disregarded for all modes. 

In other words , it is accepted that the location of the lift 
fans with respect to the craft C.G. could not affect signi- 
ficantly the roll and pitch stability. Also, air stagnation 
force in the z -direction is assumed to be negligible. 

3. During dynamic response of the craft, only air 
leakage with constant leakage area at the stern seal is 
assumed. No other air escape is considered such as possible 
leakage from sidewalls or bow seal. 

4. Plenum pressure and air mass thermodynamics are 
represented by ideal diabatic processes with the air constant 
specific heat ratio y = 1.4. 

5. Effective plenum roof area is developed as a function 
of craft weight. 

6. Planing forces are concentrated at the centroid of 
the keel longitudinal cross section through the fluid medium. 

7. Pitch and roll damping moments due to added mass are 
included in the craft dynamics , but their contribution to 
heave motion is neglected. 
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8. The thrust forces and all hydrodynamic resistive 
forces (drag) in the z-direction cancel each other, since 

at constant speed operation, these forces are in equilibrium. 
However they might contribute to the moment equations such 
as pitch moment , but considering their effective moment arms 
with respect to the C.G. are almost equal to the thrust moment 
arm, so that their contribution in the moment equations is 
also neglected. 

9. For seal modeling, vertical resistive forces due to 
hydrodynamic response are taken into account including 
corresponding effective wetted seal area. 

10. Simplified craft geometry is used for plen\im, 
buoyancy and seal modeling. 

B. COORDINATE SYSTEM AND SIGN CONVENTIONS 

1. The cartesian coordinate system is used with its 
origin located on the water line corresponding to the craft 
center of gravity. 

2. The x-axis lies on the water line with the positive 
longitudinal displacement being measured forward. 

3. The y-axis lies in the vertical cross section plane 
with the lateral dislacement being measured positive to 
starboard. 

4. The z-axis lies in the vertical cross section plane 
with the positive vertical displacement being measured 
downward. 
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5. z represents heave C.G. displacement measured from 
the origin (at water line) to the center of gravity. 

6. z is one of the craft dimensional values measured 

s 

from the keel to the C.G. 

7. Pitch is the rotational motion about the y-axis, 
angular pitch displacement is represented by 0 , positive 
pitch angle is measured upward (boat noses up) . 

8. Roll is the rotational motion about the x-axis, 

(|) angle represents angular roll displacement, positive roll 
angle is measured to starboard (boat heels to starboard) . 

9. Similarly the yaw is the rotational motion about the 
z-axis, positive yaw is measured in the clockwise direction. 

10. Positive roll, pitch and yaw moments are to be 
defined as the moments in the positive roll, pitch and 
yaw directions correspondingly. 

11. As a reference zero pitch and roll angle provide the 
x-y plane parallel to the water surface. 

The craft coordinate system is shown in figure 1. 

C. SIMPLIFIED CRAFT MODELING AND GEOMETRY 

1. General Simplified Geometry and Dimensions of the 
XR-3 Test Craft 

In figure 1 the sideview simplified geometry of the 
XR-3 test craft is shown. Also in figure 3 the top view is 
illustrated. The variables and dimensions in these figures 
are defined as follows; 



14 



C.G 
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Figure 1. Craft coordinate system and simplified geometry 
(side view) 




Figure 2. Bow seal geometry 




Figure 3. Simplified geometry (top view) 
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= Draft of the craft at center of gravity 

C = Plenum pressure center (the centroid of the 
^ effective plenum area) 

= The longitudinal distance of the sections' 
centroids from the C.G. (forward and aft 
section centroids respectively) 

* Bow and stern seals front face distances 
at water level 

= Centroid distance of the planing force 
effective area 

^l'^2 ~ Forward and aft sections sidewall submerged 

region length. 

All craft dimensional values are taken from the XR-3 
six degree-of-freedom model scale drawings as indicated in 
[2] and [3]. 

2 . Sidewalls 

As reported in [2] , the sidewall modeling is developed 
from the SIDEWALL subroutine of the 6 DOF computer model. 

Like craft dimensions, all sidewall dimensions and their 
related geometric configurations are based on scale model 
drawings. For simplicity of calculation, the sidewalls are 
divided into two sections called forward and aft. During 
dynamic response of the sidewalls, these sections are con- 
sidered individually with their own centroids, so related 
calculations can be done easily. Especially, the sidewall 
modeling is developed in order to determine buoyancy forces 
and moments. By means of this modeling, the average sidewall 
width, average draft at section centroid and effective 
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plenum area are also determined and used extensively in the 
rest of the calculations. 

The sidewall vertical cross sections at the section 
centroid forward and aft are sketched in figures 4(a) and 4(b) 
respectively. From these figures, the average forward side- 
wall width at the buoyancy center is found to be 



^f 2 tan 3* 



(II-l) 



where 

= Draft at forward buoyancy center 

tan 3^ = Tangent of the sidewall deadrise angle ' 

at the forward buoyancy center 

= Sidewall flat surface width at the kell 
(forward section) 



Therefore submerged sidewall volume of the forward 
section is given as follows 



Vi = ClI-2) 

Z^ can be formulated as 

Z^ = “ ^5 tan 9 

then the small angle approximation tan 9=0 gives 
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Forward section 
Figure 4 (a) 




Aft section 
Figure 4 (b) 
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(II-3) 



- ^5 ® 



As defined before, draft at C.G. is 



Z 



d 



z + z 



s 



(II-4) 



Similar expressions can be derived for the aft 



section , 



d , , 

2 tan 3o 



(II-5) 



where 



sidewall flat surface width at the keel 
(aft section) 



V, 



'B I 
a 2 a 



(II-6) 



+ *-6 ® 



Note that I are functions of draft and pitch 

X Si 

angle . 

3. Plenum Modeling 

In [2] plenum modeling is presented including plenum 
pressure variation due to changing of draft. In the develop- 
ment of a pitch, heave and roll linear model, the simplified 
plenum geometry and dynamics which were derived by Barnes 
are considered and used in our plenum dynamics representation. 

The variables and dimensions for plenum modeling, 
as seen in figure 5, are found as 
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(II-8) 



^8 " ^9 ^ h ^ 



A. = b Z 
^ w p 



^9 ^8 " 



(II-9) 



where 

Ap = Length of pleniim at the water line 
Aj^ = Effective pleniim area 

b = Lateral distance from one sidewall to 
another in the pleniim region. 

^ t ^oi Zj, Z^ and h are all constants related to 
9 8 d o 

craft dimensions. Only the effective plenum area varies 
when the craft draft varies, that means variation of the plemom 
pressure centroid is forward or backward. 

4 . Bow and Aft Seals 

The seals are an involved additional dynamic system 
with their own degree of freedom relative to the craft, and 
they carry little load compared with the other components 
of the craft. However in the simplified model, the seal 
motions relative to the craft are neglected. 

The assumed seal geometry is shown in figure 2, 
developed in [2] . Using an appropriate geometric configura- 
tion, the seal wetted area is derived as a function of craft 
draft and pitch angle. 
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From [2] , the width of the wetted surface for the 



bow seal is ; 



- I- tan 9 

y = Q 3 

^f sin 31“ 



(II-IO) 



where 



sin 31“ = sin of the angle between water surface 

and the bow seal front face. 

The numerator of this ratio represents draft at the 
bow seal front face. 

A similar expression for the stern seal is; 

+£. tan 9 
_ d 4 

^a sin 32“ 



Multiplying x, and x by b , bow and stern seals 

^ A VV 

wetted areas are obtained as follows; 



= b 



w 



X, 



= b 



w 



X 



5. Rudder and Propulsion Hydrodynamics 

The rudder hydrodynamics is calculated from [4] , 
inserting rudder geometry dimensional values into the fin 
equations. Actually we are interested in the lateral force 
(sway) on the rudders for a given rudder order deflection. 
This dynamic equation is given in [4] as follows; 



22 



(Il-ll(a) ) 



r 



j d 

= p u A^d + t:^) C [6^ - (1 +r^) -2-j 



r r 



hr u 



where 



= Rudder platform area (for one rudder only) 



= draft at the stern = z + z„ - x 9 

s r 

= Vertical depth to bottom tip of the rudder 

= Rudder lift coefficient 
2 it R 

= 5 — V— T- (for fully wetted flow regime) 

^a ^ 

= The rudder aspect ratio 

= Rudder deflection (the sign convention of the 
rudder is to be positive for port turn) 



V. 



V + X r - z p 
r s ^ 



X = longitudinal coordinate of rudder geometric 
center (always negative value) 

V = lateral speed (ft/sec) 

r = yaw rate (rad/sec) 

In order to simplify the rudder dynamic equation, 
the following coefficients are defined as 



d 

r 



p u^ A^(l + jp) 
r 
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Then equation (Il-ll(a)) becomes 



Y = C. (6 - C — ) (Il-ll(b)) 

r t r e u 

Let us define to be the vertical coordinate of 
the rudder from the C.G. to its geometric center, the rudder 
roll and yaw moments can be calculated from 

Kj. = - (II-ll(c)) 

^r " '^r ^r (.11-11 (,d)) 

During turn maneuvers, the propulsion thrust vector 
is also deflected according to the given rudder order. There- 
fore the lateral component of the thrust contributes yaw and 
roll moments as well as a sway mode, but this diturbance must 
be considered as a step input to the linear modes with constant 
magnitude , therefore the total propulsion sway force is 

Y = 2 sin 6 ClI-12) 

pr t r 

where is the propulsion thrust corresponding to craft speed. 

Similarly the roll and yaw moments generated by the 
propulsion system are 




Y z 
pr p 



ClI-13 (a) ) 
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(11-13 (b) ) 




Y 

pr 



X 



P 



where x and z are the longitudinal and vertical locations 
of the propeller respectively. 

D. DYNAMIC FORCES FOR FIVE DEGREE-OF-FREEDOM MODES 
1. Seal Forces 

Considering the seal configuration and dynamics , the 
pressure difference across the wetted portion of the seal is 



P 



b 



- P 

a 



where 

— 2 
Pj^ = Plenum gauge pressure (Ib/ft ) 

P = Atmospheric pressure 

Since the seal is in equilibrium, the vertical hydro- 
dynamic force on the seal is 



'sf 



“ ^b Af 



= - P, 



i, - Jl,0 

^w ^ sin 31”^ 



(11-14) 



For the stern seal hydrodynamics, a similar expression 
can be found as 



sa 



- (P. + P.)* 

b 1 d sin 32°' 



(11-15) 
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where 



P. = Additional pressure supplied into stern air 
^ bags in order to reduce leakage. 



2 . Buoyant Forces 

The forward section buoyant force is determined from 
equation (II-2) , it is 



Z 



bf 



2 p g 



here p is the density of the seawater and g is the gravitational 
constant. 

Finally, the forward and aft section forces are' 



Zbf - 2 p g(^ g 



— + fa j) 



2ba ' - 2 0 



..2^ H- t) ) £ P. 

2 tan 6. a^ a 2 



(11-16) 

(11-17) 



3 . Planing Forces 

In addition to the pitch angle 9, the sidewall flat 
area at the keel determines the planing effect of the craft. 
When the craft pitches up, the frontal area (through the fluid 
medivim) causes a lift-drag force in the z-direction. 

The planing force, from Reidel's work 13], is; 



1 2 

Z„ = --s-irpAu sin 9 
pX, 2 ^ 
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where 



A 



flat surface area at the keel 



u 



speed of the craft 



Using the small angle approximation and including the 



loss coefficient suggested by Barnes in [2] , the planing 

force is redetermined as 



Z 



pi 



2 

C , TT p A u 6 
pi 



ClI-18) 



A similar expression is derived for sidewall dynamic 
in [4] , but in our modelling equation (11-18) is used. 

4 . Pitch Damping Forces 

The sidewalls of the surface effect ship influence 
moderately heave, pitch and roll modes, so in [4] the hydro- 
dynamic and hydrostatic forces and moments acting on the side 
walls of uniform triangular cross section are derived with 
particular application to vertical plane motion and lateral 
motion. 



neglecting some terms in the equations provided from [4] , 
the vertical force acting on only one sidewall is found as; 



Considering constant speed operation of the craft 



Z 



sw 



A.. X u q 



s 



(.11-19) 



1 
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where 



A ^2 = vertical added mass calculated at the stern 

= (it/8) P 

B = value of sidewall waterline at the sidewall 
“ stern 

X = stern distance from C.G. (negative value) 
s 

q = pitch rate (rad/sec) . 



For the pitch damping moment, the same equation is 
used by Barnes in [2] . 

5. Roll Damping Forces 

The XR-3 test craft is dynamically similar to 100-B 
surface effect ship. Therefore we can use dynamic equations 
calculated in [4] for the XR-3 test craft. The heave-roll 
damping force for one sidewall can be written as; 



In the above equation, is half of the craft width 



and it must be used (+) for starboard sidewall, (-) for port 
sidewall. 

Similarly the lateral damping force due to sidewall 
hydrodynamics is 



Z 




( 11 - 20 ) 



Y 



sw 



1 



u(v + x^r - z p) 

^ S S 



ClI-21) 
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here ^22 lateral added mass calculated at the stern, 

p is the roll rate. ^22 calculated from [4] , that is 

A 22 = tt/ 5 P d^2 

6. Plenum Pressure Lift Force 

From [2] and [4] , the plenum pressure lift force is 
determined assuming uniform pressure distribution. This is 
formulated as 






from equations (II-8) and (II-9) follow 



_ ^ 

Z = - fa Pu(^'o 

p w b 8 V / 



( 11 - 22 ) 



E. DYNAMIC PLENUM AIR MASS, PRESSURE EQUATIONS 

Assuming constant aft seal leakage area, the leakage flow 
rate given by Gerba and Thaler for the XR-3 heave 

only model in [1] , then we have 



2 P 



^out 



= C A„ { 



n £ p . 



b) 1/2 



(11-24) 



where 

C = Nozzle flow coefficinet 
n 

A^ = Leakage area 
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Atmospheric air density 



The linearized air flow into the plenum is pre- 

sented in [5] according to the fan properties and calculated 
in [1] as follows. 



°in = - ’'p ''b> 

where n is the number of fans supplying the air to the plenum, 
Qq is the steady state air flow rate corresponding to steady 
gauge pressure, kp is (1/slope) of the vs. curve at 

the equilibrium point. 

From equations (11-23) and (11-24), the plenum air flow 
rate can be written as 







^a^^in ^out^ 



(11-25) 



The equation (II-25) reveals thatif the fan flow rate 
(Qin) equals to the leakage flow rate r then the plenum 

air mass (Mj^) is constant. 

The absolute plenxam pressure is calculated from 



M, Y 

p = P ( ^— ) 



(11-26) 



In the above equation, is defined as plenum volume 
and its simplified form is calculated from [IJ , 
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(11-27) 



V, 



where is to be the empty plenum volume. 

F. NET HEAVE FORCE 

The resultant heave force is the summation of the all 
external dynamic forces, it is 



Z 



Z ^ + 
sf 



Z + Z, _ + Z, + 
sa bf ba 



Z , + Z + 
pi sw 



Z 

P 



where 



Z = Z + Z 

sw SWj^ SW2 

Z also includes a pair of heave buoyancy forces due to 

w Vv 

roll, but they are in opposite directions. Therefore these 
stabilizing forces are only taken into account in the roll 
motion. 

G. PITCH AND ROLL MOMENTS 
1. Seal Moments 
a. Pitch Mode 

The seal force is modeled by Barnes in [2] acting 
on the center of the seal wetted length. The effective pitch 
moment arms are ; 



= ^3 “ ^t /2 forward seal) (11-28) 

£. = 5,. + X /2 (for aft seal) (11-29) 
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From this follows the forward and aft seal moments 



M 



sf 



'sf 




(11-30) 



M 



sa 



sa sa 



(11-31) 



where Z ^ and Z have already been defined. 

S H S3. 

b. Roll Mode 

The flexible seals can contribute little to roll 
stability. Actually the seal configurations and hydrodynamics 
are relatively more complex. In [4], the seal roll moments 
are calculated by combining the roll contribution of all seal 
elements. 

Considering the seal response, the submerged 
region of the seal can be calculated and its roll moment con- 
tribution can be formulated in the following equations; 

Ksf = - P g x^ tan (p (11-30) 

Ksa = - P g tan (11-31) 

For this modeling, the roll moment arm is considered 
as (2/3 b^) , that is the centroid of the triangle cross section 
of the submerged portion. 

2 . Buoyant Moments 
a. Pitch Mode 

The pitch buoyancy moments are calculated from 
equation (11-16) and (11-17) assuming constant moment arms 



Jig and Jl 



The pitch buoyancy moments are 



«bf = - ^b£ S 

“ba = ^ba ^5 



b. Roll Mode 

As discussed in the heave mode, the effect of the 
sidewall buoyancy due to the roll is to be a stabilizing moment. 
From figures 4 (a) and 4 (b) the buoyancy volume change can 
be approximated as 



faf + b 



Finally, roll moment is 



«bf = - 2 P g AVj 



(11-35) 



Similarly, the aft section roll moment is 
Kj^a - - 2 p g AV^ 

3. Planing Pitch Moment 

The planing effect due to the sidewall flat area at 
the keel contributes only pitch motion as follows; 



M , = Z T Jl_ 

pi pi 7 
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4 . Pitch Damping Moment 

This is calculated from equation (11-19) as 



M 



sw 



= - 2 Z 



sw. 



- 2 u X ' 

33 s 



In the above equation, the factor "2" appears because 
there are two sidewall sections. 

5 . Roll Damping Moment 

Combining equation (11-20) and (11-21) , the roll 
damping component is 



K 



sw 



2Z b - Y z+K 

SW 2 o sWj^ s wave 



(11-37) 



where K represents the additional roll damping term 

because of vertical wave generation in the roll. 

From [3] and [4] , is 



K 



wave 



= - 32/tt b 



B P 
c 



(11-38) 



Here, is the sidewall vertical added mass given as 

2 2 

B = ir/8 p £ Z, Cotan 8 . 
c a 

6. Pleniim Pitch Moment 

From equation (11-22) and defining to be the 

plenum pressure centroid, the plenum pitch moment is given 
as 
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M 



Z X 
P cp 



A, (P, - P ) X 

b b a cp 



where X is calculated as a function of craft draft in [2] 
cp 

as follows; 



X = X + (5, H.)/2h 

cp c 9 d 



X = 5-o/2 + X 

c 8 c.g. 



7. Plenum Roll Moment 

The destabilizing static plenum roll moment is 
generated when the craft rolls to one side. Since the effec- 
tive plenum pressure centroid shifts laterally to the other 
side tending to increase rolling of the craft. The plenum 
roll moment is 



K 






- P ) (z 
a' s 



- V* 



8. Rudder and Propulsion Roll Moments 

These roll moments have already been calculated 
from equation (II-ll(c)) and (11-13 (a) ) . 



H. SWAY FORCES 

1. Cross Flow Drag (Sidewall) 

During the turn maneuver, the cross flow drag is 
generated acting on the outboard sidewall. The direction 
of this force is negative while turning to the left. 



35 



Neglecting the lateral drag component inside of the sidewall, 
the net cross flow drag can be calculated from 



D 



sw 



-ip sign(V^) ± (11-39) 



where 



V 



r 



V + X r - 



s 



- V^>P 



= cross flow drag coefficient 



Considering the outboard surface of the submerged 



sidewall sections (dead rise angle effect), the lift force 
in the z-direction can be produced due to the lateral drag 
deadrise projection. In [4] , the importance of the sidewall 
deadrise angle is mentioned, however at the hydrodynamic 
section in this reference, the lift component of the drag 
force has not been derived. 

In this modeling, the lift force (Z^) and lateral 
drag cqmponent is calculated assuming both have the same drag 
coefficients (C^ = 1.28). They are 



Y 



d 



D sin 6 
sw 



ClI-40) 



Z 



d 



D cos 3 
sw 



(11-41) 
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From equation (11-40) and (11-21) , the total hydro- 



dynamic sway force of the sidewall is 



sw 



= (D + Y ) sin B + Y 



sw 



sw. 



sw. 



(11-42) 



2 . Centrifugal Force 

The centrifugal force acting on the C.G. of the craft 
is always in the outboard direction during the turn maneuver. 
Y is the centrifugal force, given as 



Y 



c 




R, the turning radius is equal to (u/r) , combining this 
with the above equation gives 



Y = - m u r (11-43) 

c 

a negative sign is used due to the coordinate system sign 
convention. Since the right turn produces (+) yaw rate, then 
Y is calculated as (-) , so outboard sway force is generated 
for this configuration. 

3. Rudder and Propulsion Sway Forces 

They are determined from equation (Il-ll(b)) and 

( 11 - 12 ) . 

The sway forces due to the seals , air plenum and 
aerodynamic drag are neglected in this modeling. 
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I. YAW MOMENTS 



1. Sidewall Yaw Moment 

The sidewall yaw dynamic equations are given in [4]. 
After some simplification, the net widewall yaw moment is 
calculated from 



N 



sw 



Ni + N2 + 



(V ■ ""sw) 



Cll-44(a) ) 



where 



= 2u (v + x^r - ZgP) A^2^ (1 - ^) 

- ir/5* p (£ , ± fa d)) ^£u (v - z p) (11-44 (b) ) 

U w O 



N- = - p C, £, (2 s, r V + s- 9 ? V (11-44 (c)) 

z a u 1 z r 



£ + X 

X = — ^ = The ratio of C.G. location of (11-44 (d)) 

sidewall bow to the craft 
length 



^ ^ 1 -3X(1 -X) 

®1 3 



(11-44 (e) ) 



S2 = — 2— (11-44 (f)) 

X„w = - 1/2 P C, A^ u^ (11-44 (g)) 

^ P 

= port wetted area surge resistive force 
= - 1/2 p C, A u^ (11-44 (h)) 

S.W j. Wg 



= starbord wetted area surge resistive 
force 
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I 




(11-45) 



[2(i^ ± b^<p) + bgU 
2 

Average sidewall flat surface width at the keel. 

In equation (11-45) ± term related to the starboard 
or port wetted surface calculations. 

2. Rudder and Propulsion Yaw Moments 

The yaw moments due to the rudder and propulsion sys- 
tem are determined from equation (Il-ll(d)) and (11-13 (b) ) 
respectively . 

The seals, air plenum and aerodynamic yaw moments are 
neglected. Actually the aerodynamic yaw moment should be taken 
into account. Since during the turn maneuver, the yaw drift 
angle is produced, thus a positive aerodynamic yaw moment 
contribution to the yaw motion has been generated. From 6 
DOF nonlinear model simulations, this contribution is determined 
at 5.2% of the total yaw moment. 

J. TOTAL PITCH MOMENT 

The summation of the all dynamic pitch moments gives the 
net pitch moment which can be written as 

M = M - + M + M- ^ + M. +M,+M +M (11-46) 

sf sa of Da pi sw p 
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K. TOTAL ROLL MOMENT 

Before calculating net roll moment, the roll moment 
equation due to the cross flow drag term must be derived. 
Including also the added mass effect to the equation (11-40) 
and (11-41) , the result is 

‘'d “ ‘°sw * <‘’0 8) (11-47) 

Combining all roll dynamic equations , the net roll 
moment is calculated as 



K = 



K , -I- K 
sf sa 



''bf ^ 



Kk, + K + 
ba sw 



K -H K -I- 
pr r 



K -I- K, 
P d 



L. TOTAL SWAY FORCE 

From equation (Il-ll(a)), (11-12), (11-42) and (11-43), 

the total sway force is obtained as 



Y 



Y + Y + Y + Y 
r pr sw c 



(11-48) 



M. TOTAL YAW MOMENT 

The summation of the sidewall, rudder and propulsion 
yaw moments gives the net turning moment about the z-axis 



N 



N + 
sw 



N -I- 
r 



N 



pr 



(,II’-49) 
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N. FIVE DEGREE-OF-FREEDOM EQUATIONS OF MOTION 

Considering equations of motion for XR-3 test craft, 
the yaw (i|;) is defined as the rotational motion of the craft 
coordinate system relative to the fixed frame and the positive 
tl> angle is measured in the clockwise direction. This sign 
convention is in agreement with the rudder order sign defini- 
tion. Since positive rudder angle corresponds to turning 
to the left, producing negative yaw moment and finally at 
the equilibrium condition, the negative yaw moment is 
generated. 

Assuming small values of pitch (6) and roll angles and 
relating heave motion to the free surface, the relationship 
between the fixed and craft coordinate system can be written. 
Eliminating the terms related to the fixed coordinate system, 
the simplified 5 DOF equations of motion are obtained as 
follows ; 



m V 



Y 



(11-50) 



m w 



Z 



(11-51) 




K 



(11-52) 




M 



(11-53) 



I 



zz 




N 



(11-54) 
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The right hand sides of the equations of motion give the 
forces and moments acting on the craft in terms of the craft 
in terms of the craft coordinate system variables. In 
equation (11-57) w represents heave acceleration and all 
other variables in these equations have already been defined. 

In order to develop a 5 DOF linear model, only the surge 
equation of motion (x-direction) is neglected. Since from 
6 DOF nonlinear simulation results, very small velocity 
reduction is observed during the turn maneuver of the XR-3 
test craft (over the speed range 20 through 30 knots) . 
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III. LINEARIZATION OF THE EQUATIONS 



A. TAYLOR SERIES EXPANSION AND TOTAL DIFFERENTIAL 

Generally, the linearization process can be performed by 
means of a Taylor series expansion using the linear part of 
the series (omitting high order terms) . The same process can 
be easily done applying a total differential to the net force 
and moment equations, since the effect of small disturbances 
to the system is to be determined. 

The total differential of a given dependent variable A 
with respect to the independent variables x, y and z is given 
as 



B. LINEARIZATION OF THE EQUATIONS OF MOTION AND AIR FLOW 
DYNAMICS 

The total differentials of the equations of motion and 
air flow dynamics have been derived as 



dA 




ClII-1) 



dw = dZ/m 



(. 111 - 2 ) 



dq = dM/l 



YY 



(III-3) 



dm. 




dv = dY/m 



(III-5) 
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4 



4 






where represents the steady state value of the plenum 

gauge pressure according to the final equilibrium condition. 

The equations (11-52) and (11-54) both involved second 
order yaw and roll acceleration terms. Therefore applying 

the total differential and solving these two equations 

• • 

for dp and dr we have 



dr = dK/I + dN/l^ (III-6) 

dp = dK/Ij^ + dN/I (III-7) 



I, I and I, are defined as 

JL 



A I I - I ‘ 

j A XX zz 



xz 



A I I - I ' 

h XX zz xz 
"xx 



A I I - I ‘ 

h XX zz xz 



zz 



(III-8 (a) ) 



(III-8 (b) ) 



(III-8 (c) ) 



Equation (III-6) and equation (III-7) contain both net 
incremental roll (dK) and yaw (dN) moments. This is simply 
a cross product inertial term effect in the equations 

of motion. 

In equation (11-59) , dPj^ is not a state variable, it can 
be considered as an output variable, therefore from equation 
(11-26) and (11-27), dP^^ must be 
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dP. 



b 



3P. 



dz 



b 




C. FIVE DEGREE-OF-FREEDOM LINEAR MODEL SENSITIVITY COEFFICIENTS 

The sensitivity coefficients related to the right hand side 
of the equations of motion are defined as the partial deriva- 
tives of the net forces and moments divided by craft mass or 
inertial terms. 

In [2] , Barnes has developed a 2 DOF (heave and pitch) 
linear XR-3 model. Therefore most of the neave and pitch 
sensitivity coefficients are directly used in order to develop 
a 5 DOF linear model. However some modification is made for 
plenum sensitivity coefficients and pitch mode. Also some 
of the steady state operating point values have been reinves- 
tigated using the 6 DOF nonlinear model. 

In all the linearized equations , the final steady state 
values of the variables have been used. 

1. Pitch and Heave Sensitivity Coefficients 



This is claculated from (11-22) where the second term 
does not appear in [2] . 

The modified form of the plenum pitch sensitivity 
coefficient due to state variable z is to be 



In [2] equation (III-40) , the sensitivity coefficient 



due to z has been redefined as 



DHPZ 




h 
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DPPBZ 



= - Zp(0) 



-1 

2h 



cp 



Pl^(O) 



fa Jl„ 
w 9 



The derivation of these modified equations is based 
upon the total differential of the pleniim lift force, it 
is 



dZ 



8Z 3Z 

— ^ dP + 

3Pk ^ 






From equations (II-9) , (11-26) and (11-27) , the total 

differentials dVj^ and dPj^ have been calculated as 



fa 

dVb = ~(A^(0) - a^(0) \ ^ ) dz 



(III-9) 



- 1 ^ 

^"'’b = Y (III-IO) 

dPj^ and dVj^ equations are also different from the equations 
derived in [2], 

The bow and aft seal pitch moments sensitivity 
coefficients (DPSFTH and DPSATH) due to 6 and plenum air 
mass equation (AMB) are rederived and their modified forms 
are used for the pitch mode in the 5 DOF linear model computer 
program. 

Considering the 5 DOF linear model state space matrix 
representation, the ^ j ' s represent "A" matrix elements. 

Therefore the a^^j values relating to the sensitivity coeffi- 
cients for heave motion are determined as 
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DZZ 



^21 

a23 = DZTH 

324 = DZDP 

a25 = DZMB 

where DZZ, DZTH, etc., have been derived in [2]. However 
their modified forms can be found in the 5 DOF linear computer 
program. 

For pitch motion, a^j coefficients are 

a.. = DTHZ 

41 

a., = DTHTH 

43 

a.. = DTHDTH 

44 

345 = DTHMB 

2 . Plenvim Air Flow Sensitivity Coefficients 

These coefficients are calculated directly from [2]. 
The a^ j ' s related to plenum air sensitivity coefficients are 



= DMBZ 



47 



DMBMB 
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3. Sway Sensitivity Coefficients 

Before deriving the sensitivity coefficients for 
sway motion, some partial derivatives related to the rudder 
and sidewall hydrodynamics must be calculated. From equation 
(Il-ll(b)) is already known, this follows the partial 
derivatives of with respect to the state variables 



V 



3Y 
1 

3p 

3Y 
1 

3v 

3Y 

1 

JT 



Ct(Ce/u) 









Clll-11) 



(III-12) 



ClII-13) 



From equation (11-39) , the sidewall cross flow drag 

terms D partials are derived in the following equations, 
s w 

Let us define first some coefficients 



di 4 I P n 



dr = i^(0) + sign(V^) (j>(0) 



The partials are 



3D 



D 



sw 



3z 



d^ sign(V^) V^^(O) 



(III-14) 
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(III-15) 



d 2 and 



where u 



where 



V (0) 

(p 1 r o 



D = 2 d^ d^ sign(V^) (0) (Z^ - il^(0)/2) (III-16) 



D = - 2 d, d V (0) 

V 1 r r 



(III-17) 



D = - 2 d, d sign(V ) V (0) X 

r 1 r ^ ' r r s 



(III-18) 



d^ are defined as 



■^2 ' 2 *33 *’ o ^ “ 



(III-19) 



d 3 = A22 u 



(III-20) 



is the craft speed assumed to be constant. 
For sway mode, the state matrix elements are 



‘96 



l(il) 



i(D, sin 6) 

m (p 



(111-21) 



^ is given from equation (III-15) . 



197 



1 ( 11 ) 

m^3p^ 



= -[R + (R +d_Z )sin 3 + d,Z ] 

xn p p 3 s 3 s 



(III-22) 
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= 

^99 

= + (D^-d3)sin B - d^] (III-23) 

a = i(H) 

^911 m'3r^ 

= ^[R^ + (D - d-x ) sin 6 - d-X -mu] 

m r r J s 3 s 

(III-24) 

4 . Roll and Yaw Sensitivity Coefficients 

Recalling equations (III-6) and (III-7) , the net 
roll and yaw moments appear in both equations of motion. In 
this case, to determine sensitivity coefficients is not easy, 
since the yaw and roll partials are involved in the dynamic 
equation with different inertial terms and coefficients. 

In order to simplify the problem, the roll and yaw 
partial derivatives are calculated individually, then combining 
these coefficients and dividing by corresponding inertial 
terms give the system sensitivity coefficients for the roll 
and yaw modes. 

a. Roll Mode 

The net roll partial derivatives are defined by 
variables. 

The partial derivative with respect to z is 



^71 



3z 



^(K . + K + 
3z sf sa 



K 



bf 



^ba> 



(III-25) 



In equation (III-26) all roll components have 
already been derived in the preceding sections . 
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From equations (11-30) and (11-31 ) , the seal 

parti als are 



'""sf _ 

3z 


(^) P g(l/sin 31») 4>(0) 


9K 

sa _ _ 

3z 


(^) P g (1/sin 32°) (|)(0) 



For buoyancy partials, let us define the following 
coefficients ; 



^1 - 


2 p g (III-26) 


^2 - 


2 p g £2 (III-27) 


^3 ■ 


b_ + b 

<— ^ 2 tin 61 


^4 - 


tan's^ f-f 


^6 - 


tan "“62 “"a 



From equations (11-35) and (11-36), the buoyancy 

partials are 



3z 


^1 ^3 

= - tin b/^5-^^3^‘^^0) 


3K, 

ba 

3z 


^2 ^4 

= - tin 



51 



Equation (11-26) follows; 



^ - -IA^(O) P^(0) +(Z^-H^(0) *(0) 

Similarly for sidewall cross flow drag partial is 

3K. 

■s— = D (6 COS 3 - z sin 3) 

dZ Z O O 

The roll partial derivatives with respect to 



0 are 



k 



73 



3K 

30 



= 36 <"'sf * "sa + "bf * 



(III-31) 



where 



3K 



sf 



30 



3K. 



sa 



30 



3K 



bf 



30 



3K 



ba 



30 



/w .3 

<23> 0 9 



, w . 3 
^2.3^ ^ 



^1 ^2 
tan 3, ^^5 



^2 ^4 
tan 3- 



(b 



(Jl3/sin 31®) 4>(0) 
g(5,^/sin 32®) 

+ b^ Jig) ())(0) 
g + b^ Jig) <() (0) 
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is 



The roll partial derivative with respect to 



II 1! 

in 


3K 

_£ 

Aj^(O) (Zg -£^(0))^ <|)(0) (III-32) 

b 



The roll partials with respect to <j> are 



^76 


3K 
3 <j) 




^(K ^ + K + K. ^ + K, + K + K,) 

3<{) ' sf sa bf ha p d 


where 


(III-33) 


3K - 
s£ 

3(J) 


= - (6^2.3) ^ P g 


3K 

sa 

3(}) 


= - (fay 2. 3)^ P g 


3 <|) 


= - b^ bg 


3K, 

ba 

3 <() 


= - ^2 ^6 


3K 

3<() 


Pt^(O) (Z^ - i^(0)) Z 
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9(j) 



D.(b^ cos 6 - sin 3) 
(p o o 



where 



The roll partials with respect to p are 



.77 



9K 

3p 



■'d ^ "r* 



(III-34) 



3K 

sw 



- A., u - 32 /tt b ^ B 

2 s 22 o c 



9K. 
c 

9p 



3Y 



SW- 



+ -5TT 
p 3p 



-) (b^ cos 3 - sin 3) 



The roll partials with respect to V are 



^79 



3V 



W <>'sw * •'d + >'r> 



(III-35) 



where 



3K. 

c 

3V 



(D^ “ *^3) ^ 



Z sin 3) + d_ X Z 
o 3 s s 



3K 

sw 



d^ Z 
3 s 
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3Kr 



C^(Vu) 



The roll partials with respect to yaw rate are 



‘711 



9K 

3r 



= * ■'r> 



(III-36) 



where 



3K 



sw 



3r 



d- X Z 
3 s s 



3K 



sw 



3r 



(D - d,X ) (6 cos 3 
r 3 s o 



Z^ sin 3) 



3K^ 

Jr 






b. Yaw Mode 

The following coefficients are defined for yaw 
rotational motion; 



ni = 2 A 22 £(1 - X) 



(III-37) 



n2 = - tt/5 p u £ 



(III-38) 



= - sign(V^) p 



(III-39) 
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(III-40) 



= 2 



=2 - ®2 



(III-41) 



= - 1/2 p Cf £ 



(III-42) 



n±j variables. 



The net yaw partial derivatives are defined by 

s . 

The yaw partials with respect to ()> are 



3N 



n 



sw 



116 



3 (j> 

3N, 



+ -^(X 
3(j) 3(|) ^ pw 



X ) b 
sw o 



(III-43) 



where 



3N^ 

3^ 



4 V(0) 4>(0) 



‘ ""sw) = - ^ ^w 



The yaw partials with respect to p are 



n 



117 



= -^{N + N ) 

3p 3p^ sw “^r^ 



(III-44) 



The N and N partial derivatives can be calculated from 
sw r ^ 

equations (11-44 (a) ) and (II-ll (d) ) . 
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The yaw moments partial derivatives with respect 



to V and r are derived similarly, the n^^j coefficients 
related to these modes are 



^119 


= II (III-45) 


"nil 


= 1^ (III-46) 



c. Roll Sensitivity Coefficients 

The roll sensitivity coefficients corresponding 
to state matrix elements a^^’s are determined from equation 



(III-7) . They are 




II 




^73 




^75 = 


^e/^k 


^76 = 


^76/^k ^ "iie/^ 


II 




^711 " 


k7ii/Ik + 



where inertial terms I cind are found from equations (III- 
8(a)) and (III-8(c)). 
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d. Yaw Sensitivity Coefficients 

For yaw motion, the sensitivity coefficients 
a^^'s have been determined from equation (III-6) . 

are listed below; 



where the inertial term is calculated from equation 
(III-8 (b)). 

D. STATE SPACE REPRESENTATION AND OUTPUT EQUATION 

The state vector of the linear 5 DOF model is shown as 

X = [dZ dW d0 dq dMj^ d(J) dp dy dV dip dr]*^ CIII~47) 

Then, the system state equation must be 

X = Ax + Bu (.111-48) 



These coefficients 



‘111 



k7i/I 



113 



k73/I 



‘115 



k,s/I 



‘116 



k76/I + 



117 



k^g/I + 



^1111 ^ 711 ^^^ 
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And the output equation is 



y = C X (III-49) 

"A" matrix elements have already been determined. However 
for "B" matrix elements, the step input to the system can 
be considered as a weight removal disturbance and for the 
roll and yaw modes, a rudder order disturbance. 

For weight removal forcing function, the input is to be 

W - W. 

‘>21 = - 1 -^ ‘‘“- S '” 

where 

= final craft weight 

= initial craft weight 

m = craft mass corresponding to the final 
weight W^. 

When the rudder order is used, sway, roll and yaw forcing 
functions are generated as follows; 

=way ‘>91 - 5“ \ * =t "r* 

roll b^, = (2 sin 6 + C 6 ) (- z /I, + x /I) 

71 t r t r r k r 

yaw b, , , = (2 X. sin 6 + C. 6 ) (x /I + z /I) 

111 t r t r r r r 
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Finally, defining state variables x^, x^, 



etc 



• / 



the state equations of the 5 DOF linear model are given 
as follows; 



^1 = 


^2 


X2 = 


^21^1 ^23^3 ^24^4 ^25^5 


^3 = 


^4 


^4 = 


^41^1 + a^3X3 + a^^x^ + a^gX^ 


^5 = 


^51^1 ^55^5 


^6 = 


X? 


II 


ayiXi + a^3X3 + a^5X5 + a^gXg + a^^x^ 



+ a^gXg + 



II 

00 


Xg 


Xg = 


^96^6 ■'■ ^97^7 ^99^9 ^911^11 


^10 


^11 
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^11 ^ 111^1 ^ 113^3 ^ 115^5 ^ 116^6 

■'■ ^ 117^7 ^ 119^9 ^ 1111^11 

All coefficients of the state equations are defined 
in the preceding sections. 

Considering outputs of the linear model, only incremental 
changes of the variables can be observed. 
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IV. LINEAR FIVE DEGREE-OF-FREEDOM 
COMPUTER PROGRAM DESCRIPTIONS 



The Digital Simulation Language CDSL/360) is used for 
the linear 5 DOF XR-3 computer program. Eleven first order 
integrator blocks are needed for five degree of freedom simu- 
lations and for all linear runs , fourth order Runge-Kutta 
with fixed interval integration method is selected. 

In addition to the craft dimensional values, initial and 
final steady state operating points are entered into the 
program for linearization purposes and plot outputs. 

The input parameters of this program are 

V = Craft speed (ft/sec) 

WF = Final craft weight (lb) 

W = Initial craft weight (lb) 

RUDANG = Rudder input (positive for port turn) 

DR = Average deadrise angle 

PBBAR = Final steady gtate value of plenum 
pressure (Ib/ft ) 

ALD = Final steady state craft draft (inches) 

PHI = Final steady state roll angle 

FYP = The propulsion system thrust corresponding 
to the craft speed (for each of propellers) 

VL = Steady state lateral speed (ft/sec) 

R = Steady state yaw rate (rad/sec) . 
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The steady state values of the variables are determined 



from 6 DOF nonlinear model simulations. 
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V. SIX DEGREE-OF-FREEDOM NONLINEAR MODEL SIMULATIONS 



In this chapter 6 DOF nonlinear model simulations are 
discussed. 

At first, some of the runs are made in order to determine 
craft initial steady state operating point corresponding to 
craft speed and thrust inputs. Since 6 DOF model weight 
transient response can be analyzed by means of entering proper 
craft initial values into the program, this is simply an 
initial condition response of the nonlinear dynamic system. 

A. WEIGHT REMOVAL TRANSIENT 

For 30 and 20 knot craft speeds, the steady state values 
of some variables are listed in Table I. 



Table I 

Steady State Values of the Weight Removal Response 



Speed 

(knt) 


Pitch 

(degree) 


^b 

(Ib/ft^) 


Draft 

(inch) 


Thrust 

(lb) 


Weight 

(lb) 


30 


0.95 


24.8 


8.6 


355.6 


6722 


30 


0.26 


24.8 


5.6 


287.8 


6050 


20 


0.45 


24.8 


6.06 


186.2 


6050 



For 30 knot craft speed, two different steady state pitch 
angles are observed corresponding to the different weights. 

In [2] for pitch transient response, 0.65° steady state pitch 
angle value has been used for both linear and nonlinear 
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simulations. On the other hand as seen in Table I, 0.26° 
and 0.95° steady state pitch values are obtained. This dis- 
agreement may be explained by noting that different dimensional 
values were used or more probably a different modified 6 
DOF computer program deck was used. 

For 6050 lb craft weight,, the steady state pitch angle 
0.26“ is in agreement with the value determined by Menzel in 
figure 6 in [6] . 

During the weight removal response, simulations of 6 DOF, 
two different stern seal base leakage areas are used. In [1] 
a smaller leakage area is assumed for the linear heave only 
model as = 0.384. However = 3.84 also has been simulated 
and it is understood that the smaller leakage area is more 
reasonable considering other steady state variables in the 
air flow dynamic equations (11-23) and (11-24) . Since at 
the steady state condition must be satisfied, that 

permits only smaller A^ for a given and steady state 
values. Therefore A^ = 0.384 is preferred and used for 
nonlinear simulations. 

B . TURNING MANEUVER 

When the 3 DOF pitch, heave and roll linear model was 
developed, the 6 DOF nonlinear model simulations were used 
in order to determine all steady state roll moments during 
the turn maneuver. In addition to the roll moments, yaw 
and sway outputs were also investigated and the results are 
tabulated in Tables II, III and IV. 
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Table II 



Steady State Roll Moments 
(rudder = 5®, 30 knot speed) 



Roll Moments Positive Negative 

Stern seal FKSS - -2.2 

Bow seal FKBS - -280. 

Sidewall FKSW - -41.5 

Rudder FKRUD 131. 

Propulsion FKP 0.057 

Aerodynamic FKAED - -22.7 

Plenum ABPB*PHI* (-Z) 215.3 - 



Steady 
(rudder = 



Yaw Moments 

Stern seal FNBS 
Bow seal FNSS 
Sidewall FNSW 
Rudder FNRUD 
Propulsion FNP 
Aerodynamic FNAED 



Table III 
State Yaw Moments 
5°, 30 knot speed) 

Positive 

3.22 

0 . 

532.65 

2.78 



Negative 



0 . 

-482.2 



- 56.41 



66 



Table IV 



Steady State Sway Forces 
(rudder = 5®, 30 knot speed) 



Sway Forces 
Sidewall FYSW 



Positive 



Negative 



-351.8 



Rudder FYRUD 



47.6 



Aerodynamic FYAED 



19.07 



Centrifugal -R*U*AM 



419.24 



From Table II, it is understood that for a given (+) 
rudder input, the sign convention of the craft dynamics was 
in agreement with turning to the left. However, a negative 
sway displacement plot is observed for this maneuver. As a 
result, in the R.H.S. subroutine of the 6 DOF program, the 
sign correction should be done for y-displacement . 

The steady state sway forces acting on the XR-3 test 
craft in a turn to the left are shown in figure 6. As seen 
in this figure and from Table IV, the rudder sway force is 
in the steady state condition. For a left turn, initially 
the rudder sway force is generated in the reverse direction 
and finally a negative steady state rudder sway force is 
produced. Further investigation is done in the RUDDER sub- 
routine dynamic equations of the 6 DOF computer model, but 
no error has been determined in these equations. Also, 
different craft speeds and different effective rudder angle 
coefficients (ENDFAC) are used but everytime the same 
configuration is observed. 
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C.G. 




Figure 6. Steady state sway forces configuration 
at the left turn 



68 



In addition to rudder dynamics , the propulsion system 
of the XR-3 contributes to the roll, yaw and sway motions. 
However this contribution is not a dominant factor compared 
with the rudder effect. 

For the XR-3 test craft propulsion system configuration, 
the direction of the thrust vector is also changed as well 
as the rudder deflection. As seen from Tables II, III and 
the simulations with different rudder order, it is understood 
that the propulsion thrust component due to the rudder order 
is not calculated properly, therefore the following modifica- 
tion is made in the PROPULSION subroutine of the 6 DOF program. 



FYS = -STHSTS*CD+THSS*CD 

where FYS represents the sway force due to the starboard 
propeller, and the others are defined in the propulsion 
subroutine program. 

C. DEADRISE ANGLE EFFECT 

As reported in [4J , the deadrise angle has a major effect 
on the surface effect ship's roll behavior and therefore is 
an important design consideration. The deadrise angle effect 
on the roll motion has been investigated by Reidel in 13] . 

In this investigation, some arbitrary coefficients (PROMO 1, 

2, etc.) are used for the vertical force component of the 
sidewall hydrodynamics. From this investigation, the conclusion 
reported by Reidel was that more damping could be obtained for 
the pitch and roll motions by changing the deadrise angle. 
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The deadrise effect is reinvestigated by simulating the 
system with different rudder orders including zero disturbance, 
then the heave forces due to sidewall hydrodynamics were 
compared to each other!* As a result, no significant sidewall 
heave force change was observed. For instance, for 0° and 
5® rudder inputs, the sidewall heave forces (FZSW) are found 
to be -568.8 and -568.05 respectively. This result is not 
acceptable, since in a turn maneuver, the cross flow drag 
is generated at the outboard sidewall, so the deadrise projec- 
tion of cross flow drag must be contributed to the sidewall 
heave force. More investigation has been made in the SIDEWALL 
subroutine of the 6 DOF program and it has been determined 
that only the force component due to added mass given by 
equation (11-21) is projected as a lift force. Actually the 
cross flow drag component must be included to this projection, 
since the frontal area for cross flow drag component is out- 
board of the submerged section of the sidewall. Therefore 
it is concluded that the XR-3 test craft 6 DOF program SIDEWALL 
subroutine does not include effectively the lift force due 
to deadrise. On the other hand, the lift force and horizontal 
drag component can have different drag coefficients, therefore 
the modification of the sidewall subroutine is left for future 
studies . 
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VI. COMPARISON OF NONLINEAR AND LINEAR MODEL RESULTS 



A. C.G. ACCELERATION (see figs. 9 and 23) 

A 10% weight removal disturbance is used for both models. 

The linear and nonlinear model showed very rapid response 
with small rise time. As seen in figure 9, the 6 DOF model 
transient is relatively more damped compared with the linear 
model response. The linear model showed a small overshoot 
but it is also sufficiently damped. Small oscillatory behavior 
is observed in both figures which is due to pitch-heave 
coupling. In general the two models agree in transient response 
and steady state behavior. 

Both of the C.G. acceleration curves were normalized by a 
factor ”g"^ the gravitational constant. 

B. PITCH TRANSIENT (see figs. 10 and 24) 

As seen in these figures, both models showed an oscillatory 
transient response with the same frequency, w = 3.12 rad/sec. 
However their steady state values are different. In fact, 
the pitch moment due to plenum pressure and the effective 
plenum pressure center play a dominant role for pitch motion. 

In [2] Barnes developed plenum modeling including effective 
plenum roof area calculations and plenum pressure centroid 
shift due to draft variation. Therefore the linear model 
plenum modeling is different from the nonlinear modeling, 
so that a greater plenum moment arm is calculated in the linear 
model. 
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The steady state pitch angle values are 0.5® and 0.75° 
for nonlinear and linear models correspondingly. 

C. PLENUM PRESSURE TRANSIENT (see figs. 7 and 21) 

As seen in these figures, both models are exactly in 

agreement with each other. In each curve, the plenum pressure 

2 

suddenly drops to 21.6 (Ib/ft ) then increases up to its 
steady state value with sufficient damping. This is expected 
behavior because the plenum volume increases when 10% weight 
is removed, obviously this disturbance causes a sudden drop in 
plenum pressure. As seen before, pitch-plenum pressure coupling 
had been observed during the exponential increasing phase of 
the plenum pressure. At the fifth second of simulation time, 
nonlinear model plenum pressure is about 24.29 while the 
linear model pressure is 24.25. 

D. DRAFT TRANSIENT (see figs. 8 and 22) 

Excellent agreement was achieved between linear and non- 
linear model draft responses. As expected, the draft decreases 
exponentially when the weight is removed from the craft. 

The draft is found to be 5.7 inches for the linear model, 

5.85 inches for the nonlinear model at the fifth second of 
simulation time. 

As a result, both models exhibited the same draft 
responses. 

E. ROLL RESPONSE 

As mentioned before, at first a 3 DOF linear model (pitch- 
heave-roll) had been developed. 3 DOF linear and 6 DOF 
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nonlinear models roll behavior comparison is done by means 
of a step input roll moment disturbance to the linear model 
and a turning command to the nonlinear model. In figure 13 
6 DOF nonlinear model roll angle response is shown for a 5° 
rudder order Cleft turn), in this simulation 1.18® outboard 
steady state roll angle is generated. From the R.H.S. sub- 
routine of the 6 DOF program, all the steady state roll moment 
components are determined as seen in Table II and the sum of 
these disturbance moments applied to the linear model roll 
mode as a step input. Roll behavior is analyzed from figure 
27. As expected, the linear model showed a greater overshoot 
and more oscillatory response compared with the nonlinear 
roll transient. Since the linear model has been disturbed 
by a step roll moment corresponding to the nonlinear model 
steady state roll moments. On the other hand, its steady 
state roll angle is found to be 1.28®. The roll steady state 
difference of two models is about 8.4%. For a 3® rudder input, 
nonlinear and linear simulations are obtained in figures 11 
and 25 with the steady state values 0.62® and 0.65® respec'<- 
tively. Hence the difference is 4.8% in the steady state 
values for these simulations. 

Figures 12 and 14 show XR-3 yaw rate plots and corres- 
ponding linear simulations are shown in figures 26 and 28. 

As a result, the 3 DOF linear model and 6 DOF program 
roll simulations are in agreement with each other within 
10% difference. 
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As discussed in the previous chapter, some problems 
are observed in the 6 DOF nonlinear program. Therefore 
developing two more modes for the 3 DOF model, the 5 DOF 
linear model has been generated. By means of this model, 
roll behavior of the XR-3 in a turn, deadrise angle and 
cross product inertial term effects have been investigated 
for 20 knot craft speed with different rudder deflections. 

The roll response can be seen in figure 29 omitting the 

cross product inertial term (I„_ = -2800.) in the equations 

xz 

of motion, applying a 10° rudder order and assiaming 70 degree 
average deadrise angle of the craft. For this particular 
test, the craft starts inboard rolling with the oscillatory 
transient then goes to 0.47° steady state roll. Using the 
same input parameters and also including the cross product 
term I in the equations of motion, smaller overshoot (18%) 
and a more damped roll response has been observed in figure 
31. As seen in this figure, at the initial phase, secondary 
overshoot has been generated. The steady state roll angle is 
0.40° and the oscillation frequency w = 3.5 rad/sec., for 
this simulation. On the other hand, the corresponding 6 DOF 
nonlinear simulation is shown in figure 17 and also a roll 
rate plot is shown in figure 18. The nonlinear simulation 
! in figure 17 shows a more damped transient but greater steady 
state roll behavior (1.32°) compared with the corresponding 
linear simulation shown in figure 31. 
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Figures 30 and 32 present linear simulation curves for 
5 and 15 degree rudder inputs respectively. Both curves 
exhibit the same linear model response discussed before. 
Finally the linear and nonlinear model steady state responses 
are summarized as follows. 



6 DOF 5 DOF 



Rudder 


Nonlinear 


Linear 


Reference 


5 


0.57 


0.24 


figures 15 and 30 


10 


1.32 


0.40 


figures 17 and 31 


15 


2.2 


0.45 


figures 19 and 32 


The deadrise angle effect is investigated by using differ 


ent deadrise angle values for 10 degree rudder input and 20 


knot craft 


speed simulations. 


The results 


are tabulated in 


Table V. 










Table 


V 






Sidewall Deadrise 


Angle Effect 




(rudder = 10®, 20 


knot speed) 








Lateral 




Deadrise 


Roll angle 


Speed 


Yaw rate 


angle (degree) (degree) 


(ft/sec) 


(rad/sec) 


45 


-1.01 


4.18 


-0.064 


50 


-0.84 


4.04 


-0.069 


55 


-0.62 


3.90 


-0.076 


60 


-0.30 


3.83 


-0.080 


65 


0.02 


3.73 


-0.085 


70 


0.40 


3.60 


-0.100 
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For the linear model, the lateral velocity and yaw rate 
plots are shown in figures 33 and 34. Both of them showed 
sufficiently damped transient response. From figure 33, 
the lateral speed steady state value is found 3.6 (ft/sec) 
and as expected for turning to the left, the negative yaw 
rate is generated as -0.1 (rad/sec). 
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X-SCRLE=1 . OOE+OO UNITS INCH. 
r-SCFILE = 5. OOE + 00 UNITS INCH. 



XR-3 30 KNOTS SPEED 
PLOT IS PLENUM PRESSURE 



Figure 7. 6 DOF Plenum pressure transient 
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X-SCflLE=l . OOE+00 UNITS INCH. 
T-SCflLE=2. OOE+00 UNITS INCH. 



XR-3 30 KNOTS SPEED 
PLOT IS Z DISPLACEMENT 



Figure 8. 6 DOF draft transient 
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X-SCflLE=l . OOE+00 UNITS INCH. 
T-SCnLE=2. OOE-02 UNITS INCH. 

XR-3 30 KNOTS SPEED 
PLOT IS C. G. flCCELERRTION 

Figure 9. 6 DOF C.G. Acceleration 
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000 002 004 006 008 010 012 




X-SCRLE=1 . OOE+00 UNITS INCH. 
T-SCflLE=2. OOE-01 UNITS INCH. 

XR-3 30 KNOTS SPEED 
PLOT IS PITCH ANGLE 

Figure 10. 6 DOF pitch angle transient 
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X-SCFILE = 5. OOE + 00 UNITS INCH. 
Y-SCRLE=2. OOE-01 UNITS INCH. 

XR-3 30 KNOTS SPEED RD=3 
PLOT IS ROLL ONGLE 

Figure 11. 6 DOF Roll angle plot, 6 = 30° 

30 knot speed ^ 
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X-SCflLE=5. OOE+00 UNITS INCH. 

T-SCFILE = 5. OOE-01 UNITS INCH. 

XR-3 30 KNOTS SPEED RD=3 
PLOT IS ROLL RRTE 

Figure 12. 6 DOF roll rate plot, 6 =3° 

30 knot speed ^ 
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000 005 010 015 020 025 030 




X-SCRLE=5. OOE+00 UNITS INCH. 
T-SCRLE=5. OOE-01 UNITS INCH. 

XR-3 30 KNOTS SPEED RD=5 
PLOT IS ROLL ANGLE 

Figure 13. 6 DOF roll angl^ plot, 6 =5“ 

30 knot speed ^ 
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X-SCRLE=5. OOE+00 UNITS INCH. 
T-SCRLE=1 . OOE+00 UNITS INCH. 

XR-3 30 KNOTS SPEED RD=5 
PLOT IS ROLL ROTE 

Figure 14. 6 DOF roll rate plot, = 5“ 

30 knot speed 
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X-SCRLE=5. OOE+00 UNITS INCH. 
Y-SCRLE=1 .OOE-01 UNITS INCH. 

XR-3 20 KNOTS SPEED RD=5 
PLOT IS ROLL RNGLE 

Figure 15. 6 DOF roll angle plot, 6 =5“ 

20 knot speed ^ 
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X-SCRLE=5. OOE+00 UNITS INCH. 
Y-3CfiLE=2. OOE-01 UNITS INCH. 

XR-3 20 KNOTS SPEED RD=5 
PLOT IS ROLL ROTE 

Figure 16. 6 DOF roll rate plot, 6 =5° 

20 knot speed ^ 
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X-SCflLE=5. OOE+00 UNITS INCH. 

Y-SCflLE=5. OOE-01 UNITS INCH. 

XR-3 20 KNOTS SPEED RD=10 
PLOT IS ROLL ANGLE 

Figure 17. 6 DOF roll angle plot, 6 = 10° 

20 knot speed 
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X-SCRLE=5. OOE+00 UNITS INCH. 

Y-SCRLE=5. OOE-01 UNITS INCH. 

XR-3 20 KNOTS SPEED RD=10 
PLOT IS ROLL RRTE 

Figure 18. 60 DOF roll rate plot, 6 = 10° 

20 knot speed ^ 
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000 005 010 015 020 025 



X-5CflLE=5. OOE+00 UNITS INCH. 

T-SCflLE=5. OOE-01 UNITS INCH. 

XR-3 20 KNOTS SPEED RD=15 
PLOT IS ROLL ANGLE 

Figure 19. 6 DOF roll angle plot/ 6 = 15® 

20 knot speed ^ 
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X-SCnLE=5. OOE+00 UNITS INCH. 

T-SCflLE=l . OOE+00 UNITS INCH. 

XR-3 20 KNOTS SPEED RD=15 
PLOT IS ROLL RATE 

Figure 20. 6 DOF roll rate plot, 6 = 15° 

20 knot speed 
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PLOT IS PLENUM PRESSURE 
XR-3 30 KNOTS SPEED 
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XSCflLE= 1.00 UNITS/INCH RUN NO. 1 

YSCRLE= 5.00 UNITS/INOH PLOT NO. 1 



Figure 21. Linear 5 DOF plenim pressure transient 
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Figure 22. Linear 5 DOF draft transient 
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Figure 23. Linear 5 DOF C.G. acceleration 
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PLOT IS PITCH ANGLE 
XR-3 30 KNOTS SPEED 
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Figure 24. Linear 5 DOF pitch angle transient 
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PLOT IS ROLL ANGLE 
XR-3 30 KNOTS SPEED RD=3 




XSCflLE= 5.00 UNITS/INCH RUN NO. i 

TSCflLE= 0.20 UNITS/INCH PLOT NO. 1 

Figure 25. Linear 3 DOF roll angle plot, 6^ = 3° 

30 knot speed 
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XSCFILE= 5.00 UNITS/INCH BUN NO. 1 

TSCFILE= 1.00 UNITS/INGH PLOT NO. 2 

Figure 26. Linear 3 DOF roll rate plot, 6^ = 3° 

30 knot speed 
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Figure 27. Linear 3 DOF roll angle plot, 6^ = 5° 

30 knot speed 
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PLOT IS ROLL RfiTE 

XR-3 30 KNOTS SPEED RD=5 




XSCRLE= 5.00 UNITS/INCH RUN NO. 1 

rSCflLE= 1.50 UNITS/INCH PLOT NO. 2 

Figure 28. Linear 3 DOF roll rate plot, 6 =5° 

30 knot speed 
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XSCRLE= 5.00 UNITS/INCH RUN NO. 1 

TSCF1LE= 1.00 UNITS/INCH PLOT NO. 1 

Figure 29. Linear 5 DOF roll angle plot, 6 =10° 

20 knot speed (omitting I . terS) 

xz 
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PLOT IS ROLL ANGLE 
XR-3 20 KNOTS SPEED RD=5 
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TSCflLE= 0.07 UNITS/INCH PLOT NO. 1 

Figure 30. Linear 5 DOF roll angle plot, 6j- = 5° 

20 knot speed 
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XSCflLE= 5.00 UNITS/INCH RUN NO. 1 

YSCflLE= 0.15 UNITS/INCH PLOT NO. 1 

Figure 31. Linear 5 DOF roll angle plot, 6 =10° 

20 knot speed ^ 
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Figure 32. Linear 5 DOF roll angle plot, 6 = 15° 

20 knot speed ^ 
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Figure 33. Linear 5 DOF lateral velocity plot 
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Figure 34. Linear 5 DOF yaw rate plot 



104 



VII. CONCLUSIONS AND RECOMMENDATIONS 



In the preceding sections, the 5 DOF linear XR-3 model 
was developed, then weight removal response and roll behavior 
are compared with 6 IDOF nonlinear simulations. Also some 
of the craft sensitivity parameters are investigated for the 
XR-3 test craft. 

The results of this study are summarized and recommenda- 
tions are made as follows. 

1. The system response is investigated by both a linear 
model and a 6 DOF nonlinear model by means of sudden weight 
removal from the C.G. and by sway force disturbance in a 
turn. 

Comparison of results shows that the difference is 
within 10% for both models. In fact, assuming accuracy in 
measurements of about 10%, the linearized model gives 
acceptable results as compared with the nonlinear model. 
Obviously the linear model simulations take less computer 
CPU time. 

2. Concerning roll response during the turn maneuver, 
the product inertia term (I„_) affects transient response 
such as overshoot, damping, etc. More overshoot and 
oscillatory response have been observed without I 

X z 

However, this term is found to be -2 800. from the 6 DOF 
program. It is suggested that this number be verified. 
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3. The importance of deadrise is observed for roll 
steady state behavior. The XR-3 test craft steady state 
roll angle depends upon the average deadrise angle of the 
sidewall. Different average deadrise angles have been 
simulated by means of the linear model and it is observed 
that the craft rolls inboard or outboard according to the 
sidewall deadrise configuration. 

The 6 DOF nonlinear model SIDEWALL SUBROUTINE 
must be modified by including a cross flow drag component 
deadrise projection (lift force) . However a new drag 
coefficient is needed for this lift component. 

4. The aerodynamic effect of the craft contributes 
about 5% of the total yaw moment. 

5. The XR-3 failure mode operation (such as fan failure) 
can be simulated by using either the linear or nonlinear 
model. But, the nonlinear model is preferable for failure 
mode simulations, since it includes sidewall and stern 
seal leakage configurations. 



I 

i 
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APPENDIX A 



COMPUTER PROGRAM LISTING 
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CUTPUT FCRCESfMOMENTS, AND RESIDUALS AT FINAL OPER POINT 
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